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APPLICATION  OP  THE  ASYMPTOTIC  METHOD  OP  CERTAIN 
PROBLEMS  OP  THE  DYNAMICS  OP  VEHICLES 

V.  A.  YaroBhevBkly 


The  present  article  generalizes  and  reduces  to  a 
compact  fom  the  rule  of  the  conservation  of  the 
’’adiabatic  Invariant.  ” 

As  an  example  of  applying  the  formulas  obtained « 
we  will  examine  problems  concerning  orbital  decay 
of  a  satellite,  the  skip-glide  motion  of  a  glider, 
and  the  oscillations  of  a  vehicle  about  the  center 
of  mass  on  a  flight  path  with  a  variable  velocity. 


1.  When  Investigating  a  niunber  of  questions  concerning  the  dynam¬ 
ics  of  vehicles,  the  problem  reduces  to  finding  a  solution  of  a  system 
of  equations  In  which  the  principal  second-order  equation  describes  a 
rapid  quaslperlodlc  motion,  and  the  parameters  of  the  system  are  var¬ 
iable,  slowly  changing  magnitudes  whose  change  Is  described  by  a  sys¬ 
tem  of  second-order  equations. 

Let  us  examine  the  system  of  equations 


(1) 

(2) 
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Here  e  Is  the  parameter  of  smallness,  t  ■  et  Is  "slow"  time, 

F  and  s^  are  differentiable  functions  of  both  arguments. 

At  constant  values  of  x  and  r^  and  when  e  -  0,  the  first  equation 
describes  a  rapid  periodic  motion. 

To  obtain  the  relationships  characterizing  the  change  of  "ainpll> 
tude"  values  of  the  variable  ^  and  the  "mean"  values  of  the  variables 
r^,  we  will  use  In  this  work  a  method  analogous  to  the  method  devel¬ 
oped  by  Kuzmak  [1,  2]. 

We  will  represent  the  unknown  functions  as  functions  of  two  var¬ 
iables,  "slow"  time  T  and  the  phase  of  the  periodic  motion  q>: 

//  -  i/,  (f .  q>)  +  ey,  (T.  (p)  +  0  (e*).  r,  »  n,  (t)  +  er„  (t,  tf)  +  0  (e*). 

^  “  “»(t)  +  e«.>,  (T)  +  0(e*). 

To  avoid  the  appearance  of  secular  members,  we  will  require  that 
the  functions  be  periodic  with  respect  to  9.  Then  the  order  of  small¬ 
ness  of  the  terms  of  the  asymptotic  expansion  will  be  retained  ln»the 
Interval  t  ~  1,  l.e.,  t  -v  l/e  [l]. 

Separating  In  all  equations  the  terms  of  the  zero  and  first  order 
of  smallness,  we  derive 


'  “•  ^  y*’  + *  i^* ^  ^ ^  y»’  + 


+ ^  + 

+  ^  y-  + o(«*)  -  0- 


(5) 


Here  9P/5y'  designates  a  partial  derivative  of  the  function  P 
with  respect  to  the  argument  e  dy/dt.  • 
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0 


(*) 


Hereafter  such  par\:lal  derivatives  of  the  function  P  will  be 
designated  as  5Fo/Sy»  5Fo/5y^  and  hTo/^ry  and  the  functions  them¬ 
selves  P(t,  rjQ,  yo,  0)  and  Sj^(t,  rj^,  yo»  o>o  Syo/^q))  as  Po  and 
Integrating  the  expression  In  e  In  Eq.  (k)  with  respect  to  q> 
from  9  to  q)  +  2t  and  taking  Into  account  that  r^^^^  Is  a  function  peri¬ 
odic  with  respect  to  q>,  we  obtain 


^-^-0.  (5) 

of 

where 

(naturally,  when  integrating  here  with  respect  to  q>  the  quantities 
T  and  rj  are  considered  parameters). 

The  term  q))  Is  determined  by  the  relation; 

1  ? 

rn  (T.  \  (s..  —  s,.I  dip,  (6  ) 

o>o  »’  '  ' 

Vflt) 

where  q>o(T)  Is  some  function  which  Is  nonessential  when  calculating 
the  zero  terns  of  the  asymtotlc  expansion. 

The  term  yo  la  a  solution  of  the  "standard"  equation 

<^1(J)~  +  Plr.ri,.,j„0]^0.  (7) 

The  value  of  o©  la  conveniently  determined  from  the  hypothesis 
that  yo(q))  -  yo(9  +  Stt). 

The  change  In  time  of  the  amplitude  values  yo,  which  Is  of  partic¬ 
ular  Interest  to  us,  cannot  be  deteirolned  from  this  equation,  there¬ 
fore  we  must  attract  the  "condition  of  periodicity"  of  the  function 
yi  [1]. 
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The  equation  (linear)  for  determining  the  function  yi  has  the 
following  foztn  (see  Eq.  (3)}: 


*'  ^  dy 


(8) 


Differentiating  Eq.  (7 )  In  terms  of  9,  we  are  easily  convinced 
that  the  term  dyo/d^  satisfies  Eq.  (8)  without  the  right  part.  Then, 
multiplying  Eq.  (7)  by  yi  and  (12)  by  dyo/b^,  subtracting  the  obtained 
relations  from  one  another.  Integrating  In  texms  of  9  from  9  to 
9  +  27r  and  taking  Into  account  that  we  easily  find  the 

first  conditions  needed  for  the  periodicity  of  the  function  yi  In 
tezns  of  9: • 


»+»» 


dFt 

6/ 


dyt 


(9) 


We  will  not  examine  the  second  condition  of  periodicity  which  en¬ 
ables  us  to  detennlne  the  correction  Oj.. 

Taking  Into  account  that  y©  Is  a  solution  of  Eq.  (7)#  we  find 

that 


where  yjj^  Is  the  amplitude  value  of  the  function  y©  Indifferently,  more 
or  less,  corresponding  to  the  vanishing  of  the  derivative  by©/b9  (the 
plus  sign  Is  ascribed  to  the  Increasing  portion  of  y©  and  the  minus 
sign  to  the  decreasing  portion  of  y©). 
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It  l8  evident  that  relation  (10)  1b  equivalent  to  the  determina¬ 
tion  of  the  "Instantaneous  an5>lltude  yjj^  for  values  y  and  dy/dt  by 
means  of  "frozen"  Eq.  (;L)  (t  -  const,  Vj  -  const,  e  ■  0). 

We  will  transform  Eq.  (9)  to  the  following  form; 


W  Ip 

+  S  S 

1—1  9  ^ 


When  Integrating  with  respect  to  <p,  the  magnitudes  r  and  r^^  must  be 
considered  as  parameters,  therefore  5yo/5q>dq)  ■  dy©.  Taking  this  Into 
account  as  well  as  relation  (10),  we  will  reduce  Eq.-  (9)  to  the  fol¬ 
lowing  form: 


dr/  *'.•  \*/’ 

—  J  (-2  ^ 

»m 

^  S  dy,  +  2  dy, «  0. 


(11) 


(The  sign  p  means  that  Integration  Is  performed  for  the  complete  per¬ 
iod  of  the  change  yo). 

Integrating  by  parts  and  taking  Into  account  that  rj^^  Is  a  periodic 
function  deteimlned  by  formula  (6),  and 


we  obtain 


(12) 


Having  differentiated  the  first  term  In  Eq.  (11),  we  derive 
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(1?) 


Substituting  (12)  and  (13)  Into  (ll)j  we  obtain 


dx 


§  ^  + 


(14) 


Here  dyQ/dt(yo)  la  the  solution  of  Eq.  (1)  with  frozen  values  of 
and  X  when  e  -  0. 

Since 


d  rrfy,-! 
dt  dr,  I  dt\' 


then  the  third  Integral  In  (14)  can  be  presented  In  the  following  form: 


dyt 


djh 

dl 


(y*) 


i, I  i[-^H 


We  will  Introduce  the  following  designations  for  averageable 

I 

values: 
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_  f  SFt  dyt 

«■. 


f  d  r  <Zi/«  1  . 

=  IT  =, 

S/,  - - --^ - -  Sto  (t,  r 


iZm). 


^'■^dyt  =  D(x,  rjt,  y„) 


Is  the  "adiabatic  Invariant"  [jl* 

Then  the  concluding  result  can  be  presented  In  the  following  form: 


^1  .  D+ vi£.l.+D^-0. 

dX  I  r/t-cotui  ^  Z*  “ 

The  values  of  rj^  are  determined  from  Eq.  (5) 


dFt 


(15) 


=  s.t  (r,  !/»>)• 

dx 

Using  these  equations  we  can  determine  the  change  In  time  of  the 
amplitude  of  the  zero  term  of  the  asymptotic  expansion  for  (l.e., 

Yjjj)  and  also  the  change  In  time  of  the  zero  terms  of  the  a83ni5)totlc 
expansion  (mean  values)  for  r^^  (l.e.#  ’^lo^*  difference  between 

the  exact  solution  and  the  zero  terras  of  the  asymptotic  expansion  In 
the  Interval  t  ~  i/e  retains  an  order  of  smallness  e. 

The  first  two  terms  of  Eq.  (15)  can  be  combined  Into  a  "unique 
derivative"  of  the  Integral  D  In  time,  which  Is  distinguished  frcan  the 
real#  complete  derivative  In  that  the  values  drj^dt  ■  Sj^  are  re¬ 
placed  by  the  values  iT~.  If  the  right-hand  parts  of  Eq.  (2)  were  to 

V  O 

depend  on  ^  and  dy/dt#  then  (when  bF/by'  »  0)  we  would  obtain  the 
known  rule  of  the  conservation  of  the  adiabatic  Invariant  D  ■  const 
[5].  Equation  (15)  generalizes  this  rule. 

The  magnitude  bFo/by'  corresponds  to  the  term  bFo/^y'  averaged 
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I 


f’ 

t 


over  a  period  with  a  weight  of  (dyo/dt)“,  the  magnitudes  Sj^,  correspond 
to  the  terms  Sj^  averaged  over  a  period  with  a  weight 


It  la  easy  to  see  that  In  the  important  special  case  where 

the  magnitude  i^,  like  hFo/^y',  Is  a  magnitude  averaged  over  a  period 
with  a  weight  of  (dyo/dt)®. 

We  note  In  conclusion  that  these  results  can  easily  he  generalized 
for  the  case  where  the  term  e  dy/dt  Is  replaced  by  the  term  ef(dy/dt). 

2.  Let  us  examine,  as  the  first  example,  the  problem  of  the 
decay  of  an  elliptical  orbit  of  a  satellite  owing  to  a  small  aerody¬ 
namic  drag. 

The  Initial  system  of  equations  Is  written  In  the  following  form 

[4]: 


n. 


(16) 


where 


dtf  m  i  y  \  I 


(17) 


Is  the  magnitude  which  Is  the  Invariant  on  a  Keplerlan  trajectory  and 
In  this  case  slowly  changes,  R  Is  the  distance  from  the  vehicle  to 
the  center  of  the  planet,  Ro  Is  the  planet  radius,  go  Is  the  accelera 
tlon  of  gravity  at  the  planet  surface,  qpi  Is  the  central  angle,  p  Is 
the  density  of  the  atmosphere,  Is  the  drag  coefficient,  S  Is  the 
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characteristic  area,  m  Is  the  satellite  mass.  The  meaning  of  the 
designations  can  be  cleared  up  from  Pig.  1.  The  perlgean  point  cor¬ 
responds  to  the  value  ^ 

Let  us  apply  the  method  cited  in  part  1  to  this  system  (16),  (17). 
Taking  Into  account  that 

^  (£/..  n.  l)  -  n)*.  (18 ) 

we  will  find  the  value  of  the  adiabatic  Invariant 

-  t,)*. 


Fig.  1. 

Equation  (15)  has  the  following  form  here; 


^  dP  A 


Hence 


‘^Ip  _  d'l  ,  . . 

ri.p  “  dip 
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(t)o  is  the  zero  term  of  the  asymptotic  expansion  for  t)  or  the  average 
value  for  tj  during  the  circling  time). 

Equation  (5)  Is  written  as 

•  dtf  “  *5^  W- 


Here 


1’ 

Cjr£^  m.— 1)> 

<r«p  “  m  If  . 

)  I  dip 


\  i  dif  J  d£  1  dif 


(19) 


(20) 


where  d§|d9  are  detemlned  by  formula  (l8). 

Let  us  extract  the  explicit  expressions  for  these  derivatives  In 
two  extreme  cases. 

a)  Considering  the  eccentricity  of  the  orbit  as  small  (?  -  §  ^ 

Jr 

Sp)*  approximating  the  dependence  of  density  on  height  H  »  R  -  Ro 
by  the  function  p  =  po  exp(-XH),  and  Introducing  the  concept  of  an 
"average"  radius  of  the  orbit  r  »  we  will  transform  Eqs.  (19)  and 
(20)  to  the  following  (dimensional)  form: 

dR,  „-R.,  [X-(r  -  /?,)!}  X 


dr 

3^’ 


.  e-^ /, IX  (r  -  /?,))  {l  +  0[x  +  0  [^^])  • 


Designating  the  parameter  X(r  -  R^^),  which  In  the  Indicated  assump 
tlon  Is  proportional  to  the  eccentricity!  In  terns  of  x,  we  derive  the 
known  result  (see,  for  example,  [5]) 
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dr 1  In  (x) 

dx  X  /,  ijr)  ’ 

* 

where  I©  and  Ii  ere  Bespel's  functions  of  the  Imaginary  argument. 

At  large  values  of  x,  Io|li  -*•  !•  Consequently,  If  the  eccentric¬ 
ity  Is  comparatively  large,  the  height  of  the  perigee  decreases 
slowly  In  comparison  with  the  average  radius  of  the  orbit  r  (and  the 
height  of  the  apogee  R.),  and  the  eccentricity  rapidly  decreases. 

At  small  values  of  x,  lollx  "*  2|x,  x  ~  exp  Xr|2,  the  orbit 
changes  to  a  circular  orbit,  the  average  radius  of  the  orbit  and  the 
height  of  the  perigee  virtually  coincide  and  they  decrease  with  the 
same  velocity. 

b)  If  the  eccentricity  of  the  original  orbit  is  large,  then  the 
main  deceleration  of  the  satellite  occurs  In  the  portion  of  the  path 
adjacent  to  the  perlgean  point,  wherein  the  role  of  this  portion  of 
the  path  grows  with  an  Increase  of  X.  The  dependence  of  density  on 
height  for  a  large  range  of  altitudes,  characteristic  for  an  orbit 
with  a  large  eccentricity,  can  be  approximated  by  the  function 

If 

•  -  I*  KicfDUH 

|)  f>p('  ,  ^ 

where  X  Is  a  slowly  changing  function  X(eH).  In  the  section  R  -  Rp« 
«  R-,f  which  Is  of  Interest  to  us,  we  can  consider  that 


P  =  PpC 


(po  =  p  (/?.).  i(y?p)). 


simplifying  the  expressions  for  derivatives  (19)  and  (20)  In 
accord  with  the  assumptions  made,  we  obtain  (again  In  a  dimensional 
form) 
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X 


^ _ C.S/p„ 


X  [i  0  [-jr  +  0  [i;j/i;^7  ]  +  “  ]) 


C,Sp, 


m2y2fli*  Xp 


X 


Wp']) 


A  similar  result  Is  obtained  by  a  different  method  In  El'yasberg's 
study  [6].  It  follows  from  the  reduced  relations  that 


rf(i/n 


l.e.,  we  obtain  an  ordinary  linear  equation  which  Is  easily  Integrated. 

Fcrr  convenience  we  can  combine  the  results  obtained  for  small  and 
large  eccentricities  and  write  them  In  the  following  form: 


df 

d\f> 


C^Sr- 

m 


<29 


C  Sr* 


where  x 


r 

•(•<p)  ^  (r  -  Rp)  (at  large  values  of  x  I©  (x)  «  —  - 


Io(x)  -  Ii(x)  - 


•/Tr(2x)^/^ 


) 


Since  Xr  has  a  order  of  100-900,  we 


can  expect  that  the  ranges  of  variability  of  the  formulas  overlap. 

5.  As  the  second  example  we  will  examine  the  motion  of  a  flying 
vehicle  with  a  large  positive  quality  (the  so-called  skip  trajectory) 

[7]. 

We  will  Introduce  the  following  assumptions: 


R  <  ^(1.  1 0  !  <  1.  ')  C,  const. 
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Then  the  equation  of  motion  can  be  rewritten  In  the  following 

form: 

iPH  C^oSV*  _  ,  v« 

*  ~  ~  iSn  ’ 

Jcr-)  C,pSVi  djl 

nr  *  m  ‘‘‘  • 

Here  H  —  height,  V  —  velocity,  d  —  the  local  angle  of  Inclination 
of  the  trajectory,  a  —  a  variable  of  the  path,  C  and  C  —  coefficients 

A  J 

of  drag  and  lift  (see,  for  example,  [7]). 

Considering  that  the  velocity  changes  negligibly  during  the  os¬ 
cillation  period  with  respect  to  H  and  that  C  changes  slowly,  we  will 
apply  formulas  (5)  and  (I5)  to  this  system  of  equations. 

We  are  easily  satisfied  that  when  calculating  the  averaged  deriv¬ 
atives  dV®|ds  and  dV®  | ds ,  the  term  2gdH|ds  does  not  have  any  effect; 

It  vanished  on  averaging.  Therefore,  this  term  Is  not  taken  Into  ac¬ 
count  when  calculating  the  amplitude  and  mean  value  of  the  zero  terms 
of  the  asymptotic  expansion  for  height  and  velocity  respectively. 

Then,  If  we  eliminate  the  variable  ds,  we  derive  In  plane  (H,  Y) 
the  following  equation  of  motion: 


'-1 


dx* 


where  x  -  In  y 


We  can  apply  to  this  equation  the  method  cited  In  Kuzmak's  study 
[1],  considering  exp(2x)-l  a  slowly  changing  coefficient,  which  Is 
valid  for  large  values  of  the  quality  and  at  velocities  smaller  than 
circular  velocity. 

Without  going  Into  detailed  .calculations,  we  can  extract  the  final 
formula  which  determines  the  dependence  of  the  amplitude  values  of 
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height  on  the  mean  velocity: 


f  ’’iBltlal 

1  h 

[  1 

\  J 

,  ij 

1) 

Here  the  argument  of  the  function  h  Is  the  ratio  of  the  maximal 
amplitude  of  density  to  the  density  corresponding  to  a  trajectory 
of  quaslstatlonary  gliding 

The  function  h  Is  determined  by  the  formula 

h{u)  ■■■  ^  |A<  —  2  1- In t/z. 

(.0 

where  zi(u)  and  ta(u)  =»  u  are  the  roots  of  the  equation  u  -  z  + 

+  In  f  =  0. 

The  graph  of  the  function  h(u)  Is  shown  In  Fig.  2.  When 
u  1  h(u)  «  2  ^  (u  -  l)^f  when  u  »  1  h(u)  ■» 

Hence  It  Is  clear  that  In  "deep"  skips  (u  »  l)  when  C  =  const, 

O’ 

the  value  of  Is  practically  constant  (analogous  to  the  case  where 

the  height  of  the  perigee  Is  at  first  almost  constant  upon  decay  of  an 
elliptic  orbit  with  a  large  eccentricity),  which  coincides  with  the 
results  obtained  earlier  [7].  As  the  velocity  decreases  h  -*  0,  l.e., 
Pmax  trajectory  becomes  a  trajectory  of  quaslstatlonary 

gliding.  Here  the  minimal  (amplltudlnal)  density  Is  associated 

with  the  maximal  relation 

The  Instantaneous  oscillation  "period"  with  respect  to  the  vel¬ 
ocity  la  determined  by  the  formula 
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■  itK- 


i  ‘f. 


Ns  htnStlbn  ¥  (tij  u~^~r  shovm  In  Fig.  2.  since 

1{  nuittbfer  6f  coinpi6t<J  oscillations  of  the  vehicle  when  Cy  ■ 
i  ll^hSB  dolife  not  bjcceed 


^  .  V/ 


f  ^>L!i±*±jL 
i  17^  2.1C, 


X  arc  sin 


’'•nd 


VRtkt 


'inltUl 


»'m<t 


6ii 


.  v;^c,_^7  5^ 

^  4  Cj,  ^ 


Fig.  2. 


4e  th  ^hd  last  example  we  will  examine  a  rapid  plane  oBclllatory 
ftdlioK  bf  ft  ballistic  vehicle  In  an  Idealized  medium  with  a  constant 
AlSAftlty  ih  the  abS^inCe  of  a  gravitational  field,  tor  sln5)llclty  we 
llli  ftXblude  from  the  examination  the  damping  terms  in  the  equation  of 
i^hl  ddbillations  of  the  vehidle  about  the  center  of  mass  (although 
ihik  eattnot  be  done  in  praotioal  calculations). 

!bhe  equations  are  written  In  the  following  fotm: 

^  =  AiffjJaV*. 

dV 

~at=-k^C.  a)  V*. 
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Here  ki  -  ka  -  ot  Is  the  angle  of  attack,  nia  Is  a  deriv¬ 

ative  of  the  aerodynamic  moment  coefficient,  Is  the  coefficient  of 
drag,  S  and  I  are  the  characteristic  surface  and  length,  J  Is  the  mo¬ 
ment  of  Inertia,  m  Is  mass  of  the  vehicle. 

If  the  coefficient  did  not  depend  on  a,  the  change  In  ampli¬ 
tude  of  oscillations  would  be  determined  by  the  asymptotic  formula 

* 

»m  —  ^ 

y — 

or 


~  const. 


(21) 


We  will  apply  formulas  (3)  and  (13)  to  the  extracted  equations. 

i 

Taking  Into  account 


D  =  mynan. 


we  obtain  equations  describing  the  changes  of  the  zero  terms  of  the 
asyn^totlc  expansion  for  velocity  Vo  and  for  the  amplitude 


dV» 

dt 


^  __Cj.  (C»)rf3 

- kyi _ kyicTM. 


da 


X  /<-«* 


dt 


i£. 

dt 


5  C^wVaj;,  — a*(io 

I  -  «*  <'» 

-“n 


Then  the  equation  relating  the  amplitude  to  the  mean  velocity 
takes  the  following  form:  ' 
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2dlna«  4--^=-rfInV',-0.  (22) 

If  •*  C^'do^'j*  we  derive  formula  (21).  We  will  approximate 

the  dependence  C^{a)  by  a  function  of  the  form: 

Cx  •=  a  +  i)  I  a  j"  (a  >0,  b  >  0). 


Then 


C,(a*)  =  a  -+ 


Cjr(«m)  =•  a  -h 


Uslng  these  formulas,  we  can  easily  Integrate  Eq.  (22)  and  ob¬ 
tain  the  relation  between  V©  and  Oj^j  In  an  explicit  form. 

We  will  limit  ourselves  to  a  qualitative  Investigation  of  this 
equation. 

When  m  >  0 


In  the  limiting  case  when  a  =  0  and  m  »  we  derive  that  this 
relation  approaches  zero,  l.e. ,  =»  const.  In  spite  of  the  fact  that 

the  velocity  decreases  monotonlcally  and  the  relative  change  of  veloc¬ 
ity  during  the  oscillation  period  with  respect  to  a  can  be  as  small  as 
desired.  The  physical  sense  of  the  last  result  Is  that  the  fall-off 
of  velocity  V  mainly  occurs  only  In  the  neighborhood  of  amplitude 
a  »  Ojjj,  and  outside  this  neighborhood  the  velocity  can  be  considered 
as  a  virtually  constant  value.  Therefore  the  motion  "Inside"  each 
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half-period  Is  conservative,  and  the  frequency  change  Jumpllke  from 
half -period  to  half-period. 

An  opposite  "degenerate"  case  can  be  obtained  If  we  assume  that 
0  >  m  >  -i  so  that  »  0)  « 

Then  the  relation  under  consideration  Is  greater  than  unity, 

l.e. ,  In  this  case  when  m -►  -1,  the  amplitude  can  Inci^ease  proportion¬ 
ally  1/V  (see  [22]).  Both  degenerate  cases  are  Illustrated  In  Pig.  3. 
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